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Abstract 

In this paper the valuation problem of a European call option in presence of both stochastic 
volatility and transaction costs is considered. In the limit of small transaction costs and 
fast mean reversion, an asymptotic expression for the option price is obtained. While the 
dominant term in the expansion it is shown to be the classical Black and Scholes solution, the 
correction terms appear at 0(e 1 ^ 2 ) and 0(e). The optimal hedging strategy is then explicitly 
obtained for the Scott's model. 

1 Introduction 

The intrinsic limitations of the Black-Scholes model in describing real markets behavior are 
very well known. Among the main assumptions underlying that model the assumptions of 
constant volatility and no transaction costs are probably the most relevant. In this paper 
we are going to consider the pricing problem of a European option in a model in which 
both proportional transaction costs are taken into account and the volatility is assumed to 
evolve according to a stochastic process of the Ornstein-Uhlenbeck type. To analyze this 
situation we shall follows a utility maximization procedure, following the seminal paper of 
M. H.A.Davis, V.G.Panas and T.Zariphopoulou [2]. 
If one uses the following utility function U: 

U(x) = 1 — exp (— 7x) , 

where 7 expresses the risk aversion of the investor, one gets, as the result of this analysis, a 
non linear PDE for the expected value of the utility-maximized wealth held in the underlying 
asset of the option. 

At this point we shall make two assumptions. First, that the process driving the volatility 
is fast mean-reverting with speed 0(e~ l ). Second, that the transaction costs are very small, 
0{e- 2 ). 
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The pricing of a European option in presence of small transaction costs was considered in [7J , 
where a correction term to the Black&Scholes pricing formula was derived. This correction 
term was found to be order 2/3 in the pricing cost. Moreover in [7J was found that the optimal 
hedging strategy consisted in not transacting when the process driving the stock price is in 
a strip around the classical Black&Scholes de/ia-hedging formula and in rebalancing the 
portfolio (selling or buying stocks) to keep the process inside the strip of no transaction. The 
width of the no transaction strip was found to be order 1/3 in the transaction costs. In [7J 
the volatility was supposed to be constant. More Recently V.I.Zakamouline studied optimal 
hedging of European options with transaction costs via a utility optimization approach in 

(PI). 

The pricing of a European option with fast mean-reverting stochastic volatility was consid- 
ered in in a series of papers (see e.g. [3], [5] and [4]). In the above mentioned papers the 
authors found the pricing formula whose leading order term is the classical Black&Scholes 
formula with averaged volatility. The correction term was order the square root of the char- 
acteristic time scale of the process driving the volatility. An optimal consumption-investment 
problem has been investigated in a paper by M.Bardi, A.Cesaroni and L.Manca ([8]) where 
a rigorous asymptotic analysis is performed and where the solution is characterized in the 
limit of fast volatility dynamics. 

More recently, M.C.Mariani, I.SenGupta and P.Bezdek ([9]) proposed a numerical approxi- 
mation scheme for European option prices in stochastic volatility models including transac- 
tion costs based on a finite-difference method. The stochastic volatility dynamics considered 
there is a slight generalization of that proposed by Hull and White |10| , since they consider a 
drift coefficient which is a general (regular) deterministic function of both the time and the 
underlying asset price, while their diffusion coefficient is linear in the instantaneous volatility. 

In the present paper we propose a different approximation method for European option 
pricing in stochastic volatility models with transaction costs, based on an asymptotic analysis 
which follows the approach pioneered by J.-P.Fouque, G. Sircar and R.Papanicolaou ([3]). The 
model we consider for the stochastic volatility dynamics is of Ornstein-Uhlenbeck type with 
a constant diffusion coefficient. This model has been originally proposed by Stein and Stein 
. We provide closed- formulas for European option prices in the limit of fast volatility and 
small transaction costs. 

The plan of the paper is the following: in Section 2 we introduce the multidimensional 
stochastic process for the dynamic of a portfolio of the writer of a European option. By 
solving a stochastic control problem, we obtain the Hamilton-Jacobi-Bellman equation. In 
Section 3, the asymptotic analysis is done assuming small transaction costs and fast mean 
reverting volatility. In section 4 the price of the option is computed and in section 5 the 
numerical results are provided and some conclusions outlined. 

For the reader's convenience, in Appendix A the source term of the equation obtained 
through the asymptotic analysis at 0(e) is calculated; in Appendix B the averages with 
respect to the Ornstein-Uhlenbeck invariant measure using the model introduced by L.Scott 
(PP) for the stochastic volatility are given; finally, in Appendix C the derivatives, which 
appears in the obtained corrected pricing formula, with respect to the stock price of the 
classical Black and Scholes model are recalled and collected. 
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2 Option pricing via utility maximization 

We suppose to have the following multidimensional stochastic process: 



dB = rBdt - (1 + \)S d~L + (I - p,)SdM (2.1) 

dy = dL-dM (2.2) 

dS = S (adt + /(z)rfW) (2.3) 

dz = i{m - z)dt + (3 (pdW + yjl - p 2 dz) . (2.4) 



In the above equations B and S are the risk- free (the "Bond") and the risky asset (the 
"Stock") respectively, r is the risk-free interest rate, a is the drift rate of the stock, A and p 
are the (proportional) cost of buying and selling a stock, / is the volatility function, which we 
shall suppose to depend on the stochastic variable z, which is sometimes called the volatility 
driving process. L and M are the cumulative number of shares bought or sold, respectively, 
in the time interval considered [0, T]. We keep the notations introduced in [2] and [7J, where 
the reader can find a detailed justification for the model for transaction costs just introduced. 
In what follows we shall always suppose f(z) to be a function bounded away from 0: 

< mi < f(z) < m,2 < oo, \/z. 

The process followed by the stochastic variable z is a Ornstein-Uhlenbeck process with aver- 
age m. The parameter £ is the rate of mean reversion volatility. 

The Brownian motions W and Z are uncorrelated and p is the instantaneous correlation 
coefficient between asset price and the volatility shocks. Usually one considers p < 0, i.e. 
the two processes are anti-correlated (e.g. when the prices go down the investors tend to be 
nervous and the volatility raises). For more details see [5] and [6]. 

We will suppose to deal with trading strategies absolutely continuous with respect to time, 
so that: 



L = / Ids , M = / mds . 

Jo Jo 

Therefore the process we are dealing with can be written in the form: 

dB = [rB - (1 + A)S1 + (1 - p)Sm] dt (2.5) 
dy = (l-m)cft (2.6) 
dS = S (adt + /(z)dW) (2.7) 

dz = £(m - z)dt + /3 (pdW + p 2 dz) . (2.8) 



The final value of a portfolio of the writer of a European option with strike price K, after 
following the strategy ir is: 

$ w (T,B*(T),y*(T),S(T),z(T)) = B*(T) + J (s {T ) <K )c(y w (T) , S (T)) + 

I(S (T)>K) [c(y n (T) - 1, S (T)) + K) . (2.9) 
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On the other hand the final value of a portfolio which does not include the option is simply: 

$ 1 (r,B 7r (T),y 7r (T),S(T),z(T))=B 7r (r) + c(y 7r (T),S(T)) . (2.10) 

We can now define the following value functions: 

Vj (B) = sup E (U (T, B 7r (T), y 7r (T), S (T), z(T)))) (2.11) 
TreT 

for j = 1, iu. Notice how this value functions depend on the initial endowment B. 
Following [2] we now define: 

Bj = inf{£ : Vj(B) > 0}. 

The fair price of the option C to avoid arbitrage, i.e. the amount of money that the writer 
has to receive to accept the obligation implicit in writing the option, will therefore be: 

C = B w -B 1 . (2.12) 

For this price the investor would in fact be indifferent between the two possibilities of going 
into the market to hedge the option, or of going into the market without the option. 
We can define the following function that will be useful in the sequel: 

^ (T, B 7r (T), y 7r (T), S (T), z (T) ) = (T, B^ (T) , y 7r (T), S (T), z(T)) — B 7r (T). (2.13) 

We have to find an equation for Vj . In what follows we shall suppress the index j and denote 
Vj with V . The problem we are dealing with is a stochastic control problem, where the 
control is the trading strategy m and 1. The equation for V is the Hamilton-Jacobi-Bellman 
equation. 

The Hamilton Jacobi Bellman equation: Suppose we have the following n-dimensional 
controlled stochastic process: 

dX = b(t, X, y)dt + a(t, X, y)dW. 
Let us define the performance functional: 

jy(s,x)=E(K(T,X(T))). 

In general one can define a performance functional that depends on the whole trajectory. In 
our case we do not need this generality. 
The infinitesimal operator associated with the stochastic process is: 

L y f = d t f + hd x J + a,ijd XiXj f, 

where 

<Hj = - (<?o T )ij ■ 

If one defines: 
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V = sup {J v : y = y(x) is a Markov control} , 

then we have the following result according to the Hamilton-Jacobi-Bellman Theorem: sup y {L y V} 
0. 

In the present case this turns out to be equivalent to the following otimization problem: 
max {{dyVj - (1 + X)Sd B Vj) I - {d y Vj - (1 - p)Sd B Vj) m+ 

0<l,m<k 

dtVj + rBd B Vj + aSdsVj + £(m - z)d z Vj + 
\[f{z)] 2 S 2 dssV j + \^d zz V 3 + PfSpd Sz v\ = 0. (2.14) 

With some analysis, as in [2], one gets that the above equation is equivalent to the following 
equation: 

max{(^- - (1 + \)Sd B Vj) , - (d y V 3 - (1 - fi)Sd B Vj) , 
dtVj + rBd B Vj + aSdsVj + £(m - z)0*V$ + 

i[/(z)] 2 S 2 d 55 V} + ^/3 2 9«^ + PfSpd Sz v\ = 0. (2.15) 

We now consider the case of the exponential utility function U{x) = 1 — exp (—72;). We note, 
just in passing, that this gives for Vj the following expression: 

Vj = l- inf {E [exp (- 7 B(T)) exp (-7*;)]} , 

where Vl/j has been previously introduced. 
In the above maximization problem let us change the variables passing Vj — > Wf. 

Vj = l-exp (-^(B + Wj)), 

where 

5 = exp [— r(T — t)] . 

Notice that with the above expression for Vj the price of the option C, as given in f|2. 12 jl . 
now becomes: 

C = W 1 -W W . (2.16) 

The maximization problem for Vj is equivalent to the following minimization problem for 
Wj : 

mm{(d y W 3 -(l + \)S) , (-d y Wj + (l-fJ,)S) , 
d t Wj - rWj + aSd s Wj + f (m - z)d z Wj + 



1 



2' 



^[/(z)] 2 5 2 dssWj-l(d s Wjf d zz Wj-Ud z W. 



1 
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3 Small transaction costs and fast mean reverting volatility: 
the asymptotic analysis 

We now suppose small transaction costs and fast mean reverting volatility. Moreover we will 
assume that the transaction costs are much smaller than the rate of mean reversion. 

x 2 c 1 a V2V 

A = fl = £ , £=-, p= —=r . 

£ yf£ 

Buying and selling costs are assumed to be the same for simplicity. 

We believe that our asymptotic assumptions are consistent with a situation where a large 
investor, facing very small transaction costs, is involved. In fact, in the empirical study [3] it 
is found that e ~ .005. For large investor, typically A < .01%. 

In absence of transaction costs and with a deterministic volatility e = 0, the investor would 
continuously trade and get a perfect hedge staying at y = y* , the "B&S" hedging strategy. 
When transaction costs are present there is a strip of small thickness around y = y* where 
he does not transact. To resolve this strip we introduce the inner rescaled coordinate Y: 

y = y * + £ a Y and d y — ► £~ a dY . (3.1) 

The unknown boundaries between the no-transaction region and the buy and sell regions are 
located at: 

y = y* + e a Y + and y = y* - £ a Y~ . 

It is very important from the practical hedger point of view to determine Y + and Y~. 
We impose the following matching conditions (see e.g. [7]): 

W NT (Y = Y ± ) = W{y = y* ± continuity 

dyW NT (Y = y ± ) = £ a d y W(y = y* ± e a y ± ) continuity of the 

first derivative 

dYYW NT (Y = F ± ) = £ 2a d yy W{y = y* ± e a y ± ) smooth pasting 

boundary condition 

These boundary conditions will force, in the asymptotic analysis below, a = 1/3. Therefore 
the strip of no-transaction will have a thickness 0(e 1//3 ). 

In the buy region (1" < Y~) we have the equation: 

(dyWj - (1 + A)5) = , (3.2) 

which solves to: 

W = (1 + \)Sy + H~ (t, S, A) . (3.3) 
In the sell region (Y > Y + ) we have the equation: 

(dyWj - (1 - \)S) = , 
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which solves to: 

W = (1- X)Sy + H + {t,S,X). 
In the no-transaction region we have the equation: 

d t Wj - rWj + aSd s Wj + -(m - z)d z Wj 



(3.4) 



1 



+ -V 2 [dzzWj-^Wj) 2 



+—uV2fS P 



dszWj - ZdsWjdsWj 



0, (3.5) 



whose solution will be find in what follows. 

3.1 The solution in the no transaction region 

As we said, in the no transaction region we use the rescaled variable Y defined by (|3.ip . The 
change of variable leads to the following transformation rules for the derivatives: 



8y - 




ds 


-> ds-e-Wyldy 


d t - 


-> dt-e-Wytdy 


d z 





We write the solution in the no transaction region in the form: 
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W NT = S(y* + eV 3 Y) + U (S, t,z)+J2 e i/6 Ui(S, t, z) + e u / 6 U u (S, t, z, Y) + ... (3.6) 



i=2 



The derivative with respect to t writes as: 

n 

d t W NT = U ot + J2 £i/6u n + £U/e Pi* ~ VtUiur) + ... 



i=2 



The derivatives with respect to S writes as: 



ii 



d s W NT = y* + U 0S + e 2 / 6 (Y + u 2S ) + Y J ^U l s + e l2l& {U l2 s-y s UiAY) + ... 

i=3 

9 

d SS W NT = U oss + Yl ^Uiss + e 10/6 {Uios + (y*s) 2 U U YY) + - 



i=2 



The derivatives with respect to z writes as: 

n 

d z W NT = U 0z + ^2e l / 6 U iz + e 12 / 6 (U 12z -y* z Uu Y ) + - 



i=2 
9 



d zz W NT = U 0zz + ^izz + e 10/6 {Uiozz + (ylfUuvy) + ... 



(3.7) 



i=2 
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The derivative with respect to S and z writes as: 

9 



dszW NT = Uosz + e i/6 UiSz + e 10/6 (UioSz + vMUiayy) + ... 



i=2 



3.2 The 0(e x ) up to 0(e x / 6 ) order equations 

To simplify the notation, and following the use in [5] and [6], we define the linear operators 
Li and the non linear operator Mc- 

C U = (m- z)U z + u 2 U zz (3.8) 
dU = -vV2p ( a y^- U * + v^lfSpXJsz (3.9) 

C 2 U = U t + ^f 2 S 2 Uss-rU + rSU s (3.10) 
N C U = -v 2l 8 (U z ) 2 (3.11) 

The 0(e~ 1 ) equation is simply: 

C U +M c U = 0. (3.12) 
The above equation can be considered an ODE in z for Uo - . 

v 2 U azz + (m - z)U 0z - v 21 - (U Qz ) 2 = . (3.13) 

o 

In [6] it is proved that the only solution of an equation of this form is a U which does not 
depend on z. The conclusion we therefore draw is that: 

Uo does not depend on z . 

The 0(e~ l / & ) equations, for i = 2, • • ■ ,5 are: 

CoUi = . 

The conclusion is: 

U i = 2, • • • , 5 does not depend on z . 

3.3 The 0(1) equation 
The O(l) equation writes: 



Uot ~ r (Sy* + Uo) + aS (y* + U os ) + C U 6 



1 



2' 



^-[f{z)YS 2 Uoss-^(y* + U s) 2 = 0. (3.14) 



S 



The above equation will be analyzed in the subsection 3.5. 



S 



3.4 The 0(e 1 ^) equation 
The 0(e 1//6 ) equation is: 

C U 7 = . 

We therefore get: 



U-j does not depend on z . 



3.5 The 0(e 2/6 ) equation 
The 0(e 2 / 6 ) equation is: 



£ U 8 + U 2t -r (SY + U 2 ) + aS(Y + U 2S ) + 



1, 



+^lf(z)] 2 S 2 U 2S s - 2^ (y* + Uos) (Y + U 2S ) = 0. (3.15) 



2' 



,7 
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In the above equation there are terms that do not depend on Y, and terms linear in Y. They 
must be equal to zero separately. ^From the terms linear in Y one gets: 

y* = -Uos + ^0- (3.16) 

The above expression gives the leading order (in absence of transaction costs) optimal hedging 
strategy. One recognizes the Black&Scholes <5-hedging strategy. 
If one inserts the above expression into the O(l) equation (|3.14p . one gets: 

C U 6 + d t U + h 2 S 2 d ss U + rSd s U - rU + ~ ( ° '"^ = . (3.17) 
2 27 f 2 

The above equation, considered as an ODE for Uq, is of the form: 

C U = X . (3.18) 
In [5] it is shown that the solvability condition for the equation (]3.18p is: 

(X> = , (3.19) 
where the average (•) is taken with respect to the Ornstein-Uhlenbeck process invariant 



measure: 



<X> = -4= / X (z)e- (m - z)2/2u2 dz. (3.20) 



Therefore the solvability condition for (j3.17j) is: 



\ x( a — r )2 1 

d t U + -a 2 S 2 d SS U - rU + rSd s U = - V ' ^ , (3.21) 



where a is the effective constant volatility: 

- 2 = (f) , 
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and r is denned as: 



1 

P 



Once one imposes to (|3.2ip the appropriate final condition, which will be different for the 
investor with option liability and the investor without it, then Uq is determined. One can go 
back to equation (|3.17|) and solve it for Uq. The solution can be written in the form: 



U 6 = U^(S,z,t) + U 6 (S,t) , 
where, Ug Z \s,z,t), the part of Uq which depends on z, has the following expression: 



(3.22) 



-CT 1 



\s\f-a*)U QSS + \ 5 ~{* 
2 27 



P f 2 



on the other hand Uq is a function that does not depend on z and that will be determined 
by the 0(e) equation in the asymptotic procedure. 

(z) 

We can get a more explicit representation for Uq , that will be useful in the next subsection. 
We first define the functions <p(z) and ip(z) as the solutions of the following problems: 



Therefore the above expression for m > can be written as: 



to 



/ 



--(- 

P \P 



(3.23) 
(3.24) 



u 



to 



1 1 6 

2 27 



(3.25) 



One can now go back to equation (|3.15p . collect all terms independent of Y and get the 
following equation: 

£0^8 + C 2 U 2 = 



The above equation is a Poisson problem of the type (|3.18|) . The solvability condition is: 

(£2) U 2 = . (3.26) 

Notice that the above equation is homogeneous in U 2 - Given that the final condition, both for 
the investor with option liability and for the investor without it, is 0, one gets the following 
conclusions: 



U 2 
Uh 







Ug(S,t) is independent of z . 
Therefore the correction to the Black and Scholes value comes up to 0(e 1 ^ 2 ) order. 
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3.6 The 0(e 3 ^) equation 

Using the expression ()3.16j) for y* , the (e 3 / 6 ) equation can be written as: 



C 2 U 3 + dUe+CoUg = 0. 



(3.27) 



Notice that in the above equation appears Uq which, until now we have derived only up to 
the function Ue(S,t), to be determined by to a higher order asymptotics. However, in (|3.27p 
Uq is hit by the operator C\, which cancels U§{S,t). 

Therefore one can consider ()3.27j) as a Poisson problem for Ug, whose solvability condition 
is: 



(C 2 ) U 3 = - {dU 6 ) 



(3.S 



The above equation is a Black and Scholes equation for U 3 with final condition and with a 
source term. We now want to rewrite the source term. 
By using for Uq the expression (I3.22D the operator C\ cancels the part not depending on z, 
and taking into account the expression (|3.25p , one can express the source term in the equation 
(PL28D as: 



- {CiU 6 

1 1 6 

t;S 2 Uoss<p + o _ ( a ~ r fi ) 

2 27 



V2 



(//) {S 3 U 0S ss + 2S 2 U 0SS ) 



(a 



r)S 2 U 0SS (^)-°-(a-rf(^ 



^ _ 5 

f) 1 

Therefore U3 solves the following Black and Scholes equation: 

1 

2 



U 3t + ^ 2 S 2 U 3SS - rU 3 + rSU 3S 



y/2 L 



(f<ff) (S 3 U 0S ss + 2S 2 U 0SS ) - (« - r)S 2 U 0S s (j 



7 



-(a 



■if/ 



with zero final data. 



(3.29) 



3.7 The 0(e 2/3 ) equation 

Since Uj does not depend on z, the 0(e 2 / 3 ) equation can be written as: 



£ 2 U 4 + C U 10 + v z (y* z yU 14 YY - jfS z Y z = . 



The equation (|3.30p can be considered as an ODE in Y for U14. It writes as: 

Uuyy = AY 2 + B , 
where we have defined the following quantities: 



(3.30) 



(3.31) 



11 



A ~ 1 ^ 



D 



5 U 2 (y .)2> 



v 2 (y*z) 

The above equation (|3.3ip solves to: 

U 14 = ^Y 4 + ^BY 2 + CY + D, (3.32) 

with C and D independent of Y. Now we have to impose the matching conditions. 
Being: 

Wbuy = (1 + e 2 )Sy + H~ (S, z, t) in the outer 

buy region 

(3.33) 

W S ell= (1 - e 2 )Sy + H + (S,z,t) in the outer 

sell region 

and imposing the continuity of the gradient at the two boundaries: 

dyW NT (Y = -Y-) = e 1,z d y W BUY {y = y* - a 1 ^-) , 
dYW NT (Y = Y+) = e l ^d y W SELL {y = y* + e V3y+) , 

one gets, at the 0(e 14 / 6 ): 

dyU u (Y = -Y~) = S , (3.34) 
d Y U u (Y = Y + ) = -S . (3.35) 

Therefore, using (|3.32|) . one gets: 

- - (y-) 3 -BY' + C = S, (3.36) 
3 

- (Y + f + BY + + C = -S. (3.37) 
o 

Moreover, being W in the outer regions linear in y, one imposes the continuity of the second 
derivative as follows: 

dyy W NT (Y = ±Y ± ) = , 

i.e.: 

A(Y + ) 2 + B = 0, 

Aiy-f + B = 0. 

Prom these equations one sees that, at this order, the bandwidth about the Black and Scholes 
strategy is symmetric, i.e.: 

Y+ = ¥-=(--) (3.38) 
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Subtracting the two equations ()3.36p and ()3.37p to eliminate C, and using the above expres- 
sions for Y , one gets: 

l^ B f/ 2 A- l l 2 = S. 



After some manipulations, and using the expressions for A and B, the equation (|3.3ip leads 
to the following equation: 



£ U l0 + C 2 U 4 



2/3 



(3.39) 



One can also find an expression for the amplitude of the no-transaction region: 



Y + = Y' 



3 1 6 



2f 2 Sj 



1/3 



(3.40) 



Equation (13. 39ft is a Poisson problem of the type of equation A3. 18f> . The solvability condition 
gives an equation for U4: 



(£2) u t 



2/3 \ 



(3.41) 



Notice also that, adding the two equations (|3.36p and (|3.37|) . one gets that C = 0. Therefore: 



c/14 = — y 4 + -BY 2 + D . 

12 2 



(3.42) 



which will be useful in the subsection 3.9. 

3.8 The 0(e 5/6 ) equation 
The 0(e 5 / 6 ) equation writes as: 



C U u+ L 2 U 5 -V2^1v P mSY-^l 



-f{z) S Y + — r [ — v 



dY 2 \ d 



0, (3.43) 



This equation can be considered as an ODE for U15: 

AY + B . 



d 2 U n 



dY 2 



where we have denoted: 



- £ Q U n + £ 2 U 5 

Integrating (|3.43p with respect to Y and using the boundary conditions: 
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u 15Y (y+) = u 15Y (-y-) = o 

which are needed to ensure the continuity of the gradient, one gets: 

C U n + C 2 U 5 = . (3.44) 
The above equation is a Poisson problem for Uu, whose solvability condition reads: 

< C 2 > U 5 = . (3.45) 

This is a homogeneous Black-Scholes equation for U5. Given that the final condition is zero, 
we gets the conclusions: 



U 5 



0, 

Un(S,t) does not depend on z . 



3.9 The 0(e) equation 
Collecting the 0(e) terms one gets: 

C 2 U 6 + dU 9 + £ U 12 - lfS 2 ^(U 3S ) 2 - uV2fSp^U 3S U 6z -y* z (m- z)U 1AY 

Id 



-^f 2 YS 2 U AS + u 2 



-y* zz U 1AY - 2y* z U UYz + (y* z ) 2 U WYY - ^(U 6z ) 2 



(3.46) 
0. 



The above equation can be considered as an ODE in Y for U\§: 

d 2 U ie 



8Y 2 



AY + B, 



where we have defined: 



A 
B 



7 f 2 S 2 U AS 
S u 2 (y 



C U 12 + CiUg + C 2 U 6 - ±f 2 S 2 ^U 2 s - uV2fSp^U 3S U 6z + 



- v 2 (y* zz U lAY + 2y* z U lAYz + ^U 2 Z ) + U 1AY y* z {m - z)] / (u 2 (y*) 2 ) . 

Notice that A does not depend on Y and in B the Y-dependent terms appear only with their 
derivatives in Y. 

We integrate the equation (|3.46p from —Y~ and Y + . 
Let us use the boundary conditions: 

U WY (Y+) = U mY (-Y-) = 0. 

Moreover, being Y + = Y~ and from the expression (|3.42p it follows that: 



/_ 



U 1AY dY = . 
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Integrating (|3.46p we get: 

C U 12 + dUg + C 2 U 6 = lf 2 S 2 ^(U 3S ) 2 + uV2fSp^U 3S U 6z + v 2l -(U %z f . 

2d 

The solvability condition for U12, gives the following equation for Uq: 

<£ 2 ) U 6 = - (Al/ 9 ) + Ia 2 5 2 ](C/ 35 ) 2 + iV&pJtfo (/^f ) + ((^ 6 ¥) 2 ) • (3-47) 

The main results of this section are the following: 

1. Equation (j33Tj) for U ; 

2. Equation (|3.26|) for C/2 which led us to U2 = 

3. Equation (IH^jl for f7 3 ; 

4. Equation (EPljl for U 4 ; 

5. Equation (|3.45p for U5 which led us to U5 = 

6. Equation (ET4T)) for C/ 6 ; 

7. The expression f|3. 16|) for y* , the center of the no-transaction region. 

8. The expression (I3.40p for the boundaries of the no transaction region. 

4 The price of the option 

To calculate the price of the option we now use the equation (12.16|) . The price will have the 
same asymptotic expansion as Wj with j = l,w, namely: 

C = C + e l ' z C 2 + v^Qs + e 2/3 C 4 + e 5/6 C 5 + eC 6 . (4.1) 

Each Cj is given by: 

a = ui - ur . 

To find the appropriate final conditions for the Cj, we write the final conditions for W 1 and 
W w . They are, respectively: 

W\T) = y(T)S(T) (4.2) 

and 

W W (T) = y{T)S{T) - max(S(T) - K, 0) (4.3) 
Given the expression (|3.6p one has that the final conditions for the the Ui are the following: 

Ul(T) = for i = 0,2,3,4,5,6. (4.4) 

Uq(T) = — max(5(T) — K, 0) (4.5) 
Uf{T) = for 1 = 2,3,4,5,6. 
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4.1 The leading order price 

To calculate the leading order price we have to calculate Uq and Uq where they both satisfy 
equation ()3.2ip . Given the respective final conditions (I4.4|) and (|4.5|) one has that: 



(T-t) 
(T-t) 



5(a-r) 2 1 

27" W 

5(a-r) 2 1 



C 



BS 



2 7 f 2 

where C BS is the classical pricing formula for a European call option, i.e.: 

C BS (S,t) = SN{d 1 )-Ke' r{T - t) N(d 2 ) , 

where: 

di J^S/K )+ _ { r^ ){ T-t) 



(4.6) 
(4.7) 



a 



Vr^t 



and N(z) is the normal cumulative distribution function. ^From the above expressions for 

C (S,t) = C BS (S,t) . 



Uq one obtains: 



4.2 The 0(e 1 ^) correction 

The equation for XJ\ and Uf is (|3.26p . a homogeneous Black and Scholes equation. In both 
cases the final condition is homogeneous. Therefore XJ\ = and = and: 

C 2 (S,t) = 0. 

4.3 The 0{e 1 / 2 ) correction 

The equation for U£ and Uf is equation (|3.29j) . in both cases with homogeneous final condi- 
tion. Using respectively the expressions (|4.6p and (|4.7p in (|3.29[) . one has that: 



u. 



ur 



(T-t)^=-(a-ry(?- 
V27 \ / 



7 \ / 



ftp') (S*Ci s s + 2S 2 C^) + 



+(a - r)S 2 ( f- ) C§§ 



(4- 



(4.9) 



Therefore: 



C 3 (S,t) = -(T-t)^= 



[ftp') {S 3 d 3 s C BS + 2S 2 d ss C 



BS\ 



BS / <? 



(4.10) 



The derivatives with respect to S of C BS are explicitly calculated in the Appendix C. More- 
over, the values of {ftp') and {tp'/f) are computed in Appendix B by using the Scott's model. 
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4.4 The 0(e 2 / 3 ) correction 

The equation for U\ and Uf is equation (|3.4ip . in both cases with homogeneous final condi- 
tion. The source term for the two problems is the same: in fact y* has the same expression 
for both problems. Therefore Uf = U\ and 

C A (S,t) = . 

4.5 The 0(e 5/6 ) correct ion 

The equation for U§ and Uf is the homogeneous Black and Scholes equation (I3.45D with zero 
final condition in both cases. Then: 

C 5 (S,t)=0. 

The correction to the price without transaction costs with stochastic volatility comes up to 
0(e). 

4.6 The 0(e) correction 



To compute the 0(e) correction we have to solve the equation (|3.47p . It is a Black and 
Scholes equation with a source term. We know that Uq is decomposed in a part dependent 
on z and a part that does non depend on z, see equation (|3.22|) . The same decomposition 
holds also for C% : 



where: 
and 



a 



u, 



(z)l 



u, 



(z)w 



c 6 = ui z)l - ui 



We have already computed Ug , as given in equation (I3.25D , we can therefore calculate . 
In fact, using (|3.25[) and the expressions (|4.6p and (|4.7p . one gets: 



U, 



(z)l 



u, 



6 

(z)m 



27 



(a — r) 2 ip 



which gives: 



a 



(z) 



1 



-b C ss ip . 



(4.11) 
(4.12) 

(4.13) 



We are now left with the task of computing Cq. 

The equation for Cq can be derived using equation (|3.47p . Subtracting the two equations 
relative to Uq and Uf one gets: 



(A)C 6 



[(^}-(A^]+^fe) 2 



+ 



(4.14) 



L 



97 



U, 



(z)l 
6z 



U, 



(z)w^ 



6z 
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This equation is a Black and Scholes equation for Cq with source term. In the Appendix A 
this source term is explicitly computed and the equation (|4.14p writes as: 



<£ 2 > Cq = (T — tfA + (T-t)B + C, 



(4.15) 



where: 



A 



D 



C 



" V Va 2 [(ftp') (S 3 Cg s + 5S 2 C## 5 + 45C##) + 



4 5 
-(a — r) 



£ ) (S*CMs + 2SC§§) 



1 2 



= -» 2 P 2 S 2 (V /> - \{a - r) [</*/) (S 3 Of/ + 8S 2 Cf/ 

+ USC§i s + 4C§§) - (a - r) (fy (S 2 Cf/ + 45C## S + 2Cf#) 
S 3 (jV) [(S 3 Cf/ + 11S 2 C 5 B / + 305Cf/ + 18C## S ) (/^') 



-(a - r) (S 2 C*f + 6SCS? + 6C#&) ^ 



(a 



^\(2SC7f# + S 2 C& s ) 



// 



x (&C§§Wf)-i{pL-r)*Wf) 



'5 



(CM) 2 {^) + i£ (a - r) 2 S 2 C## <^> 



-u 2 p 2 (a — r) 



in - ; ,5 2 Cf# ( ^ ) - (S 3 C## 5 + 2S 2 C##) / - 



~y ) — [ ti " U SSS + 2b ^SS ) \~J 

-uV [(S*Cg s + 5S 3 cMs + 4^ 2 C##) <F'/> 
-(a-r)(S 3 CMs + 2S 2 CM) (G'f)] . 
Given the homogeneous final condition, the solution of (14.150 writes as: 



c 6 = (T t)3 i + (T t)2 £ + (T — t)C\ 



We have used the fact that: 



£ 2 [ 2—Dl A + S r LJ^ B + (T — t)C 



{T - t) 2 A + (T - t)B + C + 

+ ( T ~ t)3 £ 2 A + ( T ~ t)2 £ 2 B + (T - t)C 2 C 



18 



and the last three terms are zero: 



an/iBS\ an 

c 2 [ s n ^-) = S n -^C 2 C BS = 0. 



dS n J dS n 

In Appendix C the reader can find the derivatives of C BS with respect to S up to the sixth 
order. In Appendix B the averages < • > with respect to the Ornstein-Uhlenbeck process 
invariant measure are explicitly computed using the Scott's model. 



5 Numerical results 

In this section we present the main results obtained via the asymptotic method. At first we 
plot the no transaction region for different values of the volatility, ranging from 5% up to 
60%, both in the case which does not include the option, denoted by the index 1, and in the 
case which includes the option, denoted by the index w. The volatility is chosen as in the 
Scott model, f(z) = e z . In figure 1 the curves representing the Black and Scholes strategy y* 
in the absence of transaction costs and the hedging boundaries, y = y* ± esY + , are plotted 
versus S for the first problem. From the expressions f)3. 16j) and (|3.40p it follows that these 
curves are respectively given by: 



(a — r)5 
e 2z Sj 

(a — r)S 
e 2z Sj 



±£3 



6(a 



e &z S 3 -y 3 



(5.16) 



(5.17) 



The corresponding curves in the second case are plotted in figure 2 and their equations are: 



r BS , {a-r)5 

°s +^27^' ( 5 - 18 ) 



„(j<; (a — r)S , i 



6(a — r) 2 v 2 5 3 



e 6z S' 3 7 3 



(5.19) 



' s e 2z S-f 

Both in the figures 1 and 2 the strike price is K = 0.5, the risk-free interest rate is r = 0.07, 
the drift rate of the stock is a = 0.1, the risk aversion is 7 = 1, the mean volatility a = 0.2, 
the time to expiry is 0.3 and e = 1/200. 
Finally, in figure 3 it is shown the curve representing the classical Black and Scholes price of 
a European call option with the first correction obtained at 0{e^) and the second correction 
obtained at 0{e). Here the parameters are chosen as K = 100, r = 0.04, a = 0.1, 7 = 1, 
a = 0.165, the time to expiry is 3, e = 1/200. In figure 3(a) the correlation coefficient is 
p = 0, in figure 3(b) is p = —0.2. 

We want to remark that the oscillatory behavior exhibited in figure 2b), 2c), 2d) by y 
for small values of S has been already observed by A.E.Whalley and P.Wilmott [7] also in 
models with constant volatility, although this feature seems to be more pronounced in the 
present context. Moreover the thickness of the no transaction region in presence of stochastic 
volatility seems to be bigger than in model with constant volatility. The (e) 1//3 scaling was 
also a relevant feature already established in [7] which is exhibited also by the present model. 
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z) = 0.05 




(a) 



(b) 




(c) 



(d) 



Figure 1: The no-transaction region in the case which does not include the option. The dotted 
curve represents the Black and Scholes strategy y* in the absence of transaction costs, the 
other two curves represent the hedging boundaries. See the text for the choice of parameters. 
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Figure 2: The no-transaction region in the case which includes the option. The dotted curve 
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Figure 3: The dotted line represents the classical price for a European call option; the 
dashed line represents the price including the first correction at 0(e 5 ); the continuous line 
represents the price including the first correction and the second correction 0(s). (a) p = 0; 
(b) p = -0.2. 
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